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We consider electronic properties of hollandite vanadate K2V8O16, a one-dimensional zigzag-chain 
system of tig orbitals in a mixed valent state. We first calculate the Madelung energy and obtain the 
relative stability of several charge-ordering patterns to determine the most stable one that is con- 
sistent with the observed superlattice structure. We then develop the strong-coupling perturbation 
theory to derive the effective spin-orbit Hamiltonian, starting from the triply-degenerate tig orbitals 
in the V06 octahedral structure. We apply an exact-diagonalization technique on small clusters of 
this Hamiltonian and obtain the orbital-ordering pattern and spin structures in the ground state. 
We thereby discuss the electronic and magnetic properties of K2V8O16 including predictions on the 
outcome of future experimental studies. 

PACS numbers: 71.10.-w, 71.30.+h, 75.10.-b, 71.20.Be 



I. INTRODUCTION 

In the study of strongly correlated electron systems, 
vanadium oxide has been one of the central materials. In 
particular, the discovery of the phase transition associ- 
ated with the reduction of the magnetic susceptibility in 
a mixed valence compound a / -NaV205, 1,2 together with 
subsequent experimental and theoretical studies, has es- 
tablished the novel concept of the charge-ordering (CO) 
phase transition accompanied by the spin-singlet forma- 
tion. 

Recently, Isobe et al 3 reported that, in hollandite 
vanadate K2V8O16, a metal-insulator transition (MIT) 
occurs at ~ 170 K, which is accompanied by the rapid 
reduction of the magnetic susceptibility. Below the 
transition temperature, a characteristic superlattice of 
y/2a x \pla x 2c is observed^ whereby a possible CO 
phase transition accompanied by the spin-singlet forma- 
tion was proposed. 3 




FIG. 1: Schematic representation of the double string of edge- 
shared V06 octahedra in K2V8O16. 



The crystal structure of this compound belongs to a 
group of hollandite- type phases and has a VsOie frame- 
work composed of double strings of edge-shared V06 oc- 
tahedra as shown in Fig.dJ In the view of orbital physics, 
this system may be regarded as a one-dimensional (ID) 
version of LiV02 known as a possible orbital-ordering 
(OO) system of the t2 g orbitals on the 2D triangular lat- 
tice of S = 1 spins, 4,5 as in the case of a similar hollan- 
dite vanadate Bi^VsOie. 6,7 The present system K2V8O16 



however has the average valence of V 3 T5+ and thus is in 
the mixed valent state of V 3+ : V 4+ = 3d 2 : 3d 1 = 1 : 3, 
for which quite different electronic states are expected. 
Thus, the central issue in the present system is the mech- 
anism of the MIT concerning how the highly frustrated 
spin, charge, and orbital degrees of freedom at high tem- 
peratures are relaxed by lowering temperatures and what 
type of orders is realized in the ground state. 

In this paper, we first examine the charge degrees 
of freedom of this system by calculating the Madelung 
energy; we obtain several CO patterns and their rela- 
tive stability to determine the most stable one that is 
consistent with the characteristic superlattice structure 
observed in the low-temperature phase of the material. 
We then develop the strong-coupling perturbation the- 
ory starting from the triply-degenerate t2 g orbitals in the 
V06 octahedral structure and derive the effective spin- 
orbit Hamiltonian in the approximation of neglecting the 
orbital fluctuations. We apply an exact-diagonalization 
technique on small clusters to this Hamiltonian and ob- 
tain the OO pattern and spin-spin correlation functions 
in the ground state. We thereby suggest that the state of 
local singlets of two 8 = 1/2 spins coexisting with local 
high-spin clusters of S = 3/2 should be of possible rel- 
evance with the observed electronic and magnetic states 
ofK 2 V 8 Oi 6 . 

This paper is organized as follows: In Sec. II, we cal- 
culate the CO pattern within the ionic model and derive 
the effective spin-orbit Hamiltonian based on the pertur- 
bation theory. In Sec. Ill, we calculate the orbital and 
spin structures of the derived effective Hamiltonian and 
compare the results with experiment. Summary is given 
in Sec. IV. 



II. EFFECTIVE HAMILTONIAN 

Here, we derive the effective spin-orbit Hamiltonian 
based on the strong-coupling perturbation theory, where 
the unperturbed state adopted is the lowest-energy CO 
state obtained from the ionic model in the limit of van- 
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ishing hopping parameters. 



A. Madelung energy and CO patterns 

To determine the CO pattern in the ground state of 
K2V8O16, we adopt the ionic model and calculate the 
Madelung energy, 8 the electrostatic energy of an assem- 
bly of positive and negative point charges of the ions, 
which is known to be a good measure for the relative 
stability of the spatial distributions of d electrons in 
transition- metal oxides 9,10 as well as the CO patterns of 
some organic charge-transfer salts We assume various 
d electron distributions on the V ions and compare the 
Madelung energies where we use the room-temperature 
crystal structure reported for K^VsOie^ We first cal- 
culate the Madelung energy for the system with a hy- 
pothetical uniform electron distribution, where all the V 
ions are assumed to have the valence state of V 3,75+ ; we 
find the energy of —3341.386 eV/sc where sc is the unit 
cell of the observed \[2a x y/2a x 2c superlattice. We then 
calculate the Madelung energies for various CO patterns 
with V 3+ or V 4+ and find that a number of CO pat- 
terns are more stable than the energy for the uniform 
electron distribution. The most stable CO pattern that 
is consistent with the superlattice structure observed in 
the low-temperature phase of this material is found to 
have the energy of —3370.546 eV/sc. We point out that 
the obtained low-energy CO patterns commonly have the 
CO state illustrated in Fig. [2j i.e., the ID zigzag chain 
consisting of V ions always has this CO state. Therefore, 
we hereafter assume that this CO pattern is realized in 
K2V8O16. We note that this CO pattern is different from 
the one proposed by Isobe et al. (see Fig. 4 of Ref. [I), 
which has a much higher Madelung energy —3327.885 
eV/sc. Detailed discussions on the CO patterns are given 
in Ref. 




FIG. 2: Most stable CO pattern in K2V8O16 obtained from 
the Madelung energy calculations. Only a single zigzag chain 
of V ions in the material is illustrated. 



B. Perturbation theory 

Here, we develop the strong-coupling perturbation the- 
ory starting from the ground state in the strong-coupling 
limit, i.e., the CO state shown in Fig. [2l Our starting 



high-energy Hamiltonian is of the following form: 

H = Ho + H t (1) 

Ho = V ^ TliTlj + V' ^ TliTlj 
W) [ij] 

~ C taa C l/3cr' C ipv C ioccT' 
iaa' ,a^f3 

+ U ^Tli^Tliai + U' ^ Tl^Tl^ (2) 
ia i,a^f3 

H t = ~ J2 W/3(4aa C j73<x+H.C.) (3) 

(ia,j/3),a 

where c\ ao . (ci aa ) is the creation (annihilation) operator 
of an electron at site z, orbital a, and spin a =t, !• We 
define the number operators Tii aa = c\ aa Ci acn = 
Tiia] and Tii = J2 a Tii a . V and V are the intersite 

Coulomb repulsions between nearest-neighbor and next- 
nearest-neighbor sites, respectively. Jh is the Hund's rule 
coupling, and U and U' are the intra- and inter-orbital 
on-site Coulomb repulsions, respectively. We assume the 
relation U' = U — 2 Jh throughout the paper, which is 
valid in the atomic limit. ti a jp is the hopping parame- 
ter between the orbital a on site i and orbital /3 on site 
j where G {d xy ,d yz ,d zx } in the coordinate system 
shown in Fig. [TJ We retain only the direct V-V hoppings 
between the ti g orbitals because the indirect hoppings 
via the O ions are rather small A We then have the inde- 
pendent nearest-neighbor hopping parameters t a: £5, and 
t c as shown in Fig. 03 




FIG. 3: Schematic representation of the tig orbitals on the 
ID zigzag chain of K2V8O16. Two of the four lobes for each 
of the three t^ g orbitals are drawn. 



We assume that the ground state of our Hamiltonian 
Ho should have the CO pattern shown in Fig. [21 where 
there are no doubly occupied orbitals and no V 2+ or V 5+ 
ionic states since these states are highly unrealistic in this 
material. We therefore impose the condition 

U'-2V-J H > 0. (4) 

We should also note that the unperturbed (H t = 0) states 
are spin- and orbital-degenerate of the degeneracy M = 
3 n . 3A74 . 2 3iv/4 where N is the num ber of sites. This 

degeneracy is lifted by the perturbation processes. We 
develop the second-order perturbation calculation with 



respect to ti a jp assuming that U a jp is much smaller than 
U, U' , V, and V . We thereby derive the effective spin- 
orbit Hamiltonian: 



(a) 



(b) 



(c) 



H e , = g -EiM)E <Ml TT lM V 






(5) 





where and |//) (/i,// = 1, • • • , M) are the M inde- 
pendent eigenvectors of the ground state of Ho and |n) 
is the n-th excited state of Ho. E n is the corresponding 
eigenenergy of Ho (n = denotes the ground state). 

To consider the real material, we should take into ac- 
count the effect of distortion of the V06 octahedra more 
carefully because the degeneracy of the t^g orbitals can 
be lifted. To evaluate the effect of distortion, we calcu- 
late the local symmetry of the Madelung site potential 
(j) in the point charge model for the real material. Cal- 
culated result shows that the orbital d xy (see Fig. [3]) is 
much less stable for the electron to sit on than the other 



(d) 



(e) 







orbitals d yz and d zx are; i.e. 



> 



where 

we also point out that the two orbitals d yz and d zx are 
exactly degenerate due to symmetry of the lattice. We 
therefore assume that the electrons do not occupy the 
d xy orbital in the ground state as well as in the pertur- 
bation processes. A recent NMR experiment 14 seems to 
support this assumption. 

We introduce an approximation here; because the hop- 
ping parameters U a jp take the values t a th ~ t c , we 
assume % = t c = for simplicity as in Refs.landS This 

reff 



i), (b) FM-2: the ferromag- 
Jl - did® -> djd), (c) FM- 



FIG. 4: Schematic representation of five types of the bonds 
with different exchange interaction. Electrons are located in 
the shaded lobes, (a) FM-1: the ferromagnetic bond with 
the process d 2 d) — > d}d 2 — 
netic bond with the process d\cu 3 

3: the ferromagnetic bond with the process d}a 3 , ^ % ^ 3 
d\d}j, (d) AF-1: the antiferromagnetic bond with the process 
d 2 d) — > d^dfj or d\d 2 — > d 2 dj, and (e) AF-2: the antiferro- 
magnetic bond with the process djd) — > d 2 d® or d®d 2 — > djd). 
Here, we assume that the site i denotes the site of the upper 
chain of Fig. i.e., V 3+ -V 4+ -V 3+ -V 4+ chain, and the site j 
denotes the sites of the lower chain of Fig. [2] i.e., V 4+ -V 4+ - 
V 4+ chain. 



d^ dj ^ d^ dj . 



h: 



eff 



t 2 



—4Jsi ■ Sj + cl 
t 2 

4(U' - V - J H ) 4(U' -V' + Jh) 

t 2 



c = 



4(U' - V - J H ) 4(/7 / - V + J H ) ' 



i^i 



approximation means that the terms like Hfj oc t a t c and (b) The bond FM-2 obtained with the process djd) 
ij? ff oc tbt c are all neglected, retaining only the terms 
like ij? ff oc t 2 in the second-order processes. Note that 
the orbital fluctuations are completely suppressed in this 
approximation because only two orbitals connected with 
the diagonal hopping t a come out and no off- diagonal 
hopping terms appear in the effective Hamiltonian. We 
then obtain the effective spin-orbit Hamiltonian consist- 
ing of orbital-diagonal spin-subblocks with vanishing or- 
bital off-diagonal blocks. 

In the obtained effective spin-orbit Hamiltonian, we 
have five types of the bonds with different spin exchange 
interactions as shown in Fig.[H three of them (denoted as 
FM-1, 2, 3) are the bonds with ferromagnetic exchange 
interaction due to double-exchange or Hund's rule cou- 
pling mechanisms and two of them (denoted as AF-1, 2) 
are the bonds with antiferromagnetic exchange interac- 
tion due to kinetic-exchange mechanism. Defining the 
spin-1 operator on site i as Si and spin- 1/2 operator on 
site i as s^, we have the following Hamiltonian for each 
bond shown in Fig. [4j 

(a) The bond FM-1 obtained with the process d 2 d) — > 



(7) 



(c) The bond FM-3 obtained with the process d\d 



d°dj - d]d): 



Tjet 



eff 



J = 



—4Jsi ■ Sj + cl 
t 2 

4(Z7' -2V + W - J H ) 4(E/' - 2V + V + J H ) 
St 2 t 2 



(8) 



4(U' -2V + W - J H ) 4(/7' - 2V + V + J H ) ' 



(d) The bond AF-1 obtained with the process d 2 d l - 



djd 2 



d 2 ^: 



rreff 
J 



-2JSi 
t 2 

w 
w 



cl 



d?d? or d}d 2 ^d 2 d): 



( 6 ) Hf = 2JSi-Sj+ci 



(9) 



t 2 



4(V 



-2J H ) 



2{U + W-2V- V) 2{U -U' + V' + J H ) 



4(V + 2 J H ) ' 



(e) The bond AF-2 obtained with the process djd 1 
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d?d? or 



did}: 



AJsi 



cl 



(10) 



J = — c = 



t 2 



2(/7-V) 2(U-2V + V 



Here, 1 is the unit operator and we assume that the site 
i denotes the site of the upper chain of Fig. [2j i.e., V 3+ - 
V 4 +-V 3+ -V 4+ chain, and the site j denotes the sites of 
the lower chain of Fig. [2j i.e., V 4+ -V 4+ -V 4+ chain. 
Thus, we obtain the effective spin-orbit Hamiltonian 



(ij) 



(11) 



where the sum runs over all the nearest-neighbor pairs of 
sites. Note that this effective spin-orbit Hamiltonian has 
the form of block-diagonal in the spin0orbit space; i.e., 
orbital off-diagonal blocks are all zero. In other words, 
we have several 00 patterns for the iV-site systems, and 
each of them, we have the spin Hamiltonian. By diag- 
onalizing all the spin Hamiltonians, we can obtain the 
eigenstates of our effective spin-orbit Hamiltonian. In 
particular, from the lowest-energy eigenstate obtained, 
we determine the 00 pattern of the ground state in the 
parameter space. 



(a) 



0.8 



0.6 



^0.4 



0.2 



CO unstable 

■,U'- 2V-J H <0 



OO in (b) 



0.1 0.2 0.3 0.4 

VI XT 



1 3 5 7 9 11 13 15 




FIG. 5: (a) Ground- state phase diagram of the effective spin- 
orbit Hamiltonian. The CO pattern shown in Fig. [2] is un- 
stable in the shaded region, (b) Schematic representation of 
the calculated OO state in the lower-left region of the upper 
panel (a). 



III. RESULTS OF CALCULATION 

In this section, we calculate the orbital and spin struc- 
tures in the ground state of the effective spin-orbit Hamil- 
tonian Eq. ([TT]) by using the exact-diagonalization tech- 
nique on small clusters, whereby we discuss its electronic 
and magnetic properties and compare them with exper- 
iment. We assume U' = U — 2 Jh throughout the calcu- 
lations. 



A. Orbital ordering 

We here use the 16-site cluster with four 5=1 spins 
and twelve 5 = 1/2 spins (corresponding to the filling 
of 20 electrons), where the spins are coupled with the 
exchange interactions derived in Sec. II. The periodic 
boundary condition is used. Calculations are made for all 
possible OO patterns that are consistent with the two- 
fold periodicity along the c-axis observed in experiment; 
i.e., we assume the unit cell in the presence of the orbital 
ordering, which contains the 4 sites along the chain direc- 
tion (e.g., sites 1 — 4 in Fig. [5] (b)). Hereafter we assume 
the relation V = 0.6V estimated from the experimental 
interatomic distances between V ions. 

The calculated result is shown in Fig. O where we find 
that the unique orbital state illustrated in Fig.[5](b) with 
the total-spin quantum number Stot = is realized in 
the entire parameter space shown in Fig. [5] (a) unless 
the CO state is unstable. The parameter space shown in 



Fig. [5] (a) may be presumed to contain physically realis- 
tic values for vanadium oxide materials. 15 Note that this 
orbital state is stabilized by maximizing the number of 
the FM-1 bonds as shown below. 

In Fig. [6l we show calculated values of the exchange 
coupling constant J and constant c in the effective spin- 
orbit Hamiltonian Eqs. (6)- (10). We find that the value 
of c is predominantly lower for the FM-1 bond in the 
entire parameter space shown in Fig. [5] (a). Thus, we 
understand that the OO pattern shown in Fig. [5] (b) is 
stabilized by maximizing the number of this FM-1 bond. 
We should however note that, although in the small V/U' 
region the value of J for the FM-1 bond is larger than 
those of the other bonds, the values of J for the anti- 
ferromagnetic bonds are much larger than those for the 
ferromagnetic bonds in the large V/U 1 region. In partic- 
ular, a large value of J for the AF-2 bond in the large 
V/U' region is essential for the formation of the local 
spin-singlet of the two s = 1/2 spins as shown below. We 
thus find that the spin structure can be quite different 
in the different regions in Fig. [5] (a) although the unique 
OO pattern shown in Fig. [5] (b) is stabilized. 



B. Spin correlations 

We calculate the spin-spin correlation functions (Si • 
Sj), (Si • 8j), and (si • Sj) for the ground state of the 
16-site cluster of the effective spin-orbit Hamiltonian to 
discuss the spin structures of the system. The results are 
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FIG. 6: (Color online) Calculated values of the exchange cou- 
pling constant J (left panels) and the constant c (right panels) 
given in units of t^/U' . 



shown in Fig. where we in particular show (Si • Si) for 
i G s = 1/2 and (S\-Si)/2 for z E 5 = 1 in the left panels 
of Fig. [3 and (s 2 • S;) for z G s = 1/2 and (s 2 • £i)/2 for 
i G 5 = 1 in the right panels of Fig. [71 

We find in Fig. [71 that the behavior of the spin corre- 
lations is very different between two regions of the OO 
phase in the white triangle area of Fig. 5 (a): i.e., the 
regions Jh ^> V and V ^> Jh- This result suggests that, 
although the OO pattern is the same and the total spin 
$tot = in the two regions, two different spin structures 
can be realized depending on the parameter values. 

In the region Jh ^> V, we find the results that can 
roughly be described by the situation where the local 
high-spin clusters of S = 5/2 formed by four V ions 
are coupled antiferromagnetically; we find the values 
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FIG. 7: Spin- spin correlation functions calculated for the 
ground state of the 16-site cluster. The site index i is de- 
fined in Fig. [5](b); the symbols o and • indicate the sites with 
S = 1 and s = 1/2, respectively. Left and right panels rep- 
resent the spin correlations from the site 1 (S = 1) and site 
2 (s = 1/2), respectively. Upper and lower panels represent 
the spin correlations in the parameter regions Jh ^> V and 
V ^> Jh, respectively. The arrow in the lower right panel 
indicates the formation of the local spin-singlet of the two 
s = 1/2 spins. 



(S 1 -s 2 ) = 0.482, (Si-sis) = 0.435, and (S r s 16 ) = 0.498, 
which should be 0.5 if the formation of the high-spin 
cluster of S — 5/2 were complete. We also find the 
values (s 2 ■ Si)/2 = 0.241, (s 2 ■ s 15 ) = 0.194, and 
(s 2 ■ siq) = 0.234, which should be 0.25 if the forma- 
tion of the high-spin cluster of S = 5/2 were complete. 
We also find the oscillations of ±0.5 (±0.25) between two 
clusters of four spins in the upper-left (upper-right) panel 
of Fig. 7. Thus, our system in this parameter region can 
be regarded as the state of the antiferromagnetically fluc- 
tuating local high-spin clusters of S = 5/2. This state is 
illustrated schematically in Fig. [8] (a). 

In the region V ^> Jh, we find the results that can 
roughly be described by the situation where the local 
spin-singlet states are formed between two s = 1/2 spins; 
we find the value (s 2 -Ss) = —0.734, which is only slightly 
larger than the value —0.75 that is expected when the 
spin-singlet formation is complete. We also note that 
the values of (s 2 • s^) and (s 2 • Si) are very small for all 
i except i = 3, which is consistent with the formation 
of the local spin-singlet state. We should however find 
the value (Si • Si6) = 0.498, which indicates the forma- 
tion of the high-spin cluster of S = 3/2 between the two 
spins at sites 1 and 16. We also note that the oscilla- 
tion of ±0.5 between clusters of the two spins appears. 
Thus, there remains antiferromagnetic correlations be- 
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tween the high- spin clusters of S = 3/2. This correlation 
occurs because the spin-singlet formation is strong but 
not perfect. Thus, our system in this parameter region 
can be regarded as the state of the local spin-singlets of 
two 8 = 1/2 spins coexisting with the antiferromagneti- 
cally fluctuating local high-spin clusters of S = 3/2. This 
state is illustrated schematically in Fig. [51(b). 

The spin structures thus obtained are compared with 
available experimental data in the next subsection. 



(a) 



(b) 




FIG. 8: (Color online) Schematic representation of the spin 
structures, (a) The state of the antiferromagnetically fluc- 
tuating local high-spin clusters, (b) The state of the local 
spin-singlets of two s = 1/2 spins (indicated by solid circles) 
coexisting with the local high-spin clusters of S = 3/2. 



C. Comparison with experiment 

In the present theory, we have presumed that the MIT 
observed in K2V8O16 is caused by the CO of ^-electrons 
on the V ions. We suppose this to be quite a natural in- 
terpretation because of the various features observed in 
experiment; 3 e.g., the characteristic superlattice struc- 
ture and lattice distortion observed below the transition 
temperature. We have also assumed that the electrons do 
not occupy the d xy orbital (see Fig. [3|) as can be justified 
from the calculation of the Madelung site potential. This 
is also consistent with the recent NMR measurement 14 of 
the symmetry axis of this material at room temperature. 
These results first of all support the validity of our effec- 
tive spin-orbit Hamiltonian derived in Sec. II. 

Then, as for the magnetic aspects of K2V8O16, the 
rapid reduction of the magnetic susceptibility observed 
below the transition temperature seems to suggest the 
opening of the spin gap. It seems however that a small 
but finite value of the susceptibility remains finite at low 
temperatures as can be seen in Fig. 2 of Ref. 3 even after 
subtracting the Curie term coming from the presence of 
impurities. This situation seems to suggest that the for- 
mation of nonmagnetic local spin-singlets occurs below 
the transition temperature but that there still remain in- 
teracting magnetic spins. 3 The broad nonmagnetic spec- 
trum in the insulating region observed in a recent NMR 
experiment 14 might also be interesting in this respect. 



We therefore argue that these experimental situations 
are of possible relevance with the spin structure shown 
in Fig. [8] (b) where the local spin-singlets of two 5 = 1/2 
spins coexist with the antiferromagnetically interacting 
local high-spin clusters. 



IV. SUMMARY 

We have studied the electronic and magnetic proper- 
ties in hollandite vanadate K2V8O16, a possible charge 
and orbital ordering system with the mixed valent state 
of V ions with 3d 2 : 3d 1 = 1:3 and with the ti g or- 
bitals of V ions aligned on the ID zigzag chains. First, 
we have calculated the Madelung energy of the system 
and obtained the most stable CO pattern that is con- 
sistent with the superlattice structure observed in ex- 
periment. Then, by using the second-order perturbation 
theory starting from the triply-degenerate t<i g orbitals in 
the V06 octahedral structure, we have derived the ef- 
fective spin-orbit Hamiltonian. Here, we have evaluated 
the effect of distortions of the V06 octahedra from the 
local symmetry of the Madelung site potential, which 
justifies the assumption that the electrons do not occupy 
the d xy orbital. Within the approximation of neglect- 
ing the small off-diagonal hopping parameters, we have 
found that the Hamiltonian is block-diagonal with van- 
ishing orbital off-diagonal sectors. We then have used the 
numerical exact-diagonalization technique on small clus- 
ters and have obtained the orbital-ordering pattern in the 
ground state. We have also calculate the spin-spin corre- 
lation functions and have found that, depending on the 
parameter values, either the state of the antiferromag- 
netically fluctuating local high-spin clusters or the state 
of the local spin-singlets of two 5 = 1/2 spins coexisting 
with the local high-spin clusters of S = 3/2 is realized. 
By comparing these results with available experimental 
data which are quite limited at present, we suggest that 
the latter state can be in agreement with the electronic 
ground state of hollandite vanadate K2V8O16. 

Because our study presented here contains series of the- 
oretical predictions on the outcome of future experimen- 
tal studies, we hope that the present study will help one 
understand the nature of the charge, orbital, and spin 
degrees of freedom of this intriguing material. 
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